A smoother approach to scaling by suppressing monopoles and vortices 
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Suppressing monopoles and vortices by introducing large chemical potentials for them in the Wilson action for 
the SU{2) lattice gauge theory, we study the nature of the deconfinement phase transition on N% x N T lattices 
for N T —4, 5, 6 and 8 and N a — 8-16. Using finite size scaling theory, we obtain ui = y/u — 1.93 ± 0.03 for 
N T = 4, in excellent agreement with universality. The critical couplings for N T = 4, 5, 6 and 8 lattices exhibit 
large shifts towards the strong coupling region when compared with the usual Wilson action, and suggest a lot 
smoother approach to scaling. 



1. INTRODUCTION 

Universality of lattice gauge theories implies 
that actions which differ merely by irrelevant 
terms in the parlance of the renormalization 
group lead to the same continuum quantum field 
theory. While many numerical simulations have 
been performed for the Wilson action for pure 
gauge theories, other choices, some motivated by 
the desire to find a smoother continuum limit, 
have also been employed. This usually has lead to 
results in accord with universality. Investigations 
of the deconfinement phase transition for mixed 
actions, obtained by extending the Wilson action 
by addition of an adjoint coupling term, showed 
JTJQ], however, surprising challenges to the above 
notion of universality. Several finite temperature 
investigations have established the presence of a 
second order deconfinement phase transition for 
the Wilson action. Its critical exponents have 
been shown ]3| to be in very good agreement with 
those of the three dimensional Ising model, as 
conjectured by Svetitsky and Yaffe j|. This ver- 
ification of the universality conjecture strength- 
ened our analytical understanding of the decon- 
finement phase transition which, however, came 
under a shadow of doubt by the results for the 
mixed actions. For two different actions, it was 
found HU^] on a range of temporal lattice sizes 
that for large adjoint couplings a] the deconfine- 
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mcnt transition became definitely of first order 
and b] no other separate bulk transition existed. 
In fact, simulations on large symmetric lattices 
even suggested || a lack of a bulk phase tran- 
sition at an adjoint coupling where a first order 
deconfinement transition for a lattice of temporal 
size four was clearly seen. 

Recently it was shown HQ that suppression 
of some lattice artifacts such as monopoles 
and vortices do restore universality for one of the 
two actions: no first order deconfinement tran- 
sition was found in the entire coupling plane in 
that case. Here we address this question for the 
other action in the same manner and find that 
one gains an additional bonus. The approach to 
scaling seems to become smoother than that for 
the original Wilson action. 

2. ACTION AND RESULTS 

Both actions H|| have an adjoint coupling (3 a 
in addition to the usual Wilson coupling (3. The 
(3 = axis in each case describes an SO (3) model 
which has a first order phase transition. At Pa 
= oo, the theories reduce to a Z 2 lattice gauge 
theory with again a first order phase transition 
at /3 crlt w 0.44. The first order transitions extend 
into the (/?, /3a) coupling plane, meet at a triple 
point, and continue as a single line of first order 
transitions which ends at an apparently critical 
point, say, D in both cases. The proximity of D 
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to the Pa = line has commonly been held re- 
sponsible for the abrupt change from the strong 
coupling region to the scaling region for the Wil- 
son action. Following Refs. (10], one can also 
suppress Z2 monopoles and Z2 electric vortices 
for the Bhanot-Creutz action with suitable chem- 
ical potentials to restore universality in its entire 
coupling plane. Using sign(Tr F C/p) to identify 
the Z2 monopoles and vortices, one can define 
the corresponding action for such suppressions as 
below: 



Sbc,s = ^ P 1 



p 



' c 



(1 - o-c) 



(1) 



where cr c = ILeac sign(Tr F [/p) and ct ; = 
Ilpea; sign(Tr F J7p). It is clear that the above 
theory flows to the same critical fixed point in 
the continuum limit for all Pa, Am and Xe and 
has the same scaling behavior near the critical 
point. It has to be stressed though that univer- 
sality has to be tested afresh for eq. (||) to be sure 
that the above naive argument about the Am and 
Xe terms being irrelevant is correct. This is what 
we do in the following by determining a critical 
index of the deconfinement phase transition. We 
then check whether the passage to scaling is af- 
fected by studying the deconfinement transition 
as a function of the temporal lattice size. We 
chose Pa =0, Am = 1 and Xe = 5 throughout. 

We studied the deconfinement phase transition 
on N% x N T lattices by monitoring its order pa- 
rameter and the corresponding susceptibility for 
N T = 4, 5, 6 and 8 and N a = 8, 10, 12, 14, 15, 
and 16. We used the simple Metropolis algorithm 
and tuned it to have an acceptance rate ~ 40 %. 
The expectation values of the observables were 
recorded every 20 iterations to reduce the auto- 
correlations. Errors were determined by correct- 
ing for the autocorrelations and also by the jack 
knife method. The observables monitored were 
the average plaquette, P, and the absolute value 
of the average of the deconfinement order param- 
eter \L\. We also monitored the susceptibilities 



Table 1 

The values of P at which long simulations were 
performed on N% x 4 lattices, p c and the height 
of the L| -susceptibility peak, XllT*- 





P 


Pc,N a 


..max 


8 


1.37 


1.366(7) 


9.34(7) 


10 


1.344 


1.360(5) 


14.34(11) 


12 


1.331 


1.345(2) 


20.44(29) 


14 


1.34 


1.343(2) 


27.48(64) 


Table 2 








Same as 


Table 1 but for on N% 


x 5 lattices. 


N a P 


Pc,N„ 


v max 
X-\L\ 


v max 
Aprcdictcd 



10 1.545 1.570(5) 13.17(17) 

15 1.545 1.558(2) 29.57(77) 28.78(41) 

for both \L\ and P: x\l\ = N%({L 2 ) - {\L\) 2 ) and 
X p = 6N*N T ((P 2 ) - (P) 2 ). 

2.1. N T = 4 

The deconfinement phase transition on N T = 4 
lattices was studied by first making short hystere- 
sis runs on the smallest lattice to look for abrupt 
or sharp changes in the order parameter (\L\). 
In the region of its strong variation, longer runs 
were made to check whether the |L|-susceptibility 
exhibits a peak. Histogramming technique was 
used to extrapolate to nearby couplings for do- 
ing this. A fresh run was made at the x\l\ peak 
position and the process repeated until the in- 
put coupling for the run was fairly close to the 
output peak position of the susceptibility. The 
same procedure was used for the bigger lattices 
also but by starting from the p c of the smaller 
lattice. Typically 2-4 million Monte Carlo itera- 
tions were used to estimate the magnitude of the 
peak height and the peak location for each N a . 
Table 1 lists our final results for all the N a used. 
Fitting the peak heights in Table 1 to AN£, we 
obtained u> — 1.93 ± 0.03, in excellent agreement 
with the values for both the standard Wilson ac- 
tion Q and the 3-dimensional Ising model. 

2.2. N T = 5, 6 and 8 

For larger N T , we used many longer runs in the 
region of strong variation of (\L\) to obtain the 
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Table 3 

Same as Table 1 but for on iV? x 6 lattices. 





a 




..max 


v max 
^predicted 


12 


1.75 


1.735(5) 


16.34(45) 




15 


1.70 


1.702(2) 


24.39(41) 


25.13(70) 



susceptibility directly and used the histogram- 
ming only for the finer determination of the criti- 
cal coupling. Our results for (\L\) as a function of 
f3 clearly show a deconfinement phase transition 
for N T = 5, 6 and 8. This is also evident in the 
X\l\ determinations, shown in Fig. [I] for N T = 6. 
Tables 2 and 3 list the estimated maxima of X\L\ 
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the deconfinement phase transition thus verified 
on three different temporal lattice sizes, but it 
also confirms that the same physical phase tran- 
sition is being simulated on them, thus approach- 
ing the continuum limit of a — > in a progressive 
manner by keeping the transition temperature T c 
constant in physical units. 

2.3. Scaling of T c 
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Figure 2. 1/N T as a function of (3 C . The full 
lines depict the 2-loop asymptotic scaling rela- 
tion, while the dashed line denotes eq.@, nor- 
malized at N T = 8 in all cases. 



Figure 1. The |L| -susceptibility as a function of 
P on lattices with N T = 6. The continuous lines 
are extrapolations using the histogramming tech- 
nique. 

for N T = 5 and 6 for two different spatial volumes 
along with the corresponding peak locations. Us- 
ing our value for u>, and the peak height for the 
smaller spatial volume, the ^ max on the bigger 
lattice can be predicted. These predictions are 
listed in the respective tables in the last column 
and can be seen to be in very good agreement 
with the direct Monte Carlo determinations. We 
also determined lu from the peak heights in Tables 
2 and 3 and found it to be in good agreement with 
the N T — 4 value. Not only is the universality of 



Fig. | shows aT c = N~ x as a function of the 
corresponding critical (3 for both our simulations 
(squares) with suppression of monopoles and vor- 
tices and the standard Wilson action (circles). 
The latter are taken from the compilation of Ref. 

The dashed line describes a 'phenomenolog- 
ical' scaling equation, 



aT n = 



1 



4&o 



-bi/b 2 



exp 



Ab 



(2) 



which is similar to the usual scaling relation, 
shown in Fig. |^ by full lines, but with twice 
the exponent. All lines were normalized to pass 
through the N T — 8 data points. 

One sees deviations from asymptotic scaling for 
both the Wilson action and our action with sup- 
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pression of monopoles and vortices. The devia- 
tions for the same range of N T seem larger for our 
action but then one is also considerably deeper in 
the strong coupling region of the Wilson action 
where one a priori would not have even expected 
any scaling behaviour. As the agreement of our 
results with the dashed line of eq.(||) in Fig. || 
shows, scaling may hold in this region of (3 for 
the suppressed action, since the relation between 
a and (3 in this region is similar to the asymp- 
totic scaling relation, differing only in the expo- 
nent which will cancel in dimensionless ratios of 
physical quantities. It is clear that as (3 — > oo, 
the difference between the two actions must van- 
ish. The data in Fig. |^ do show such a trend 
although the limiting point is still not reached by 
N T — 8. It seems likely though that the trend 
of evenly spaced transition points for our action 
will continue and the dashed line traced by its 
transition points will merge with the Wilson ac- 
tion by N T ~ 25 or so. If this were to be so, 
a much smoother approach to continuum limit is 
to be expected after the suppression of monopoles 
and vortices. In particular, one expects that di- 
mensionless ratios of physical quantities at the 
deconfinement phase transition couplings should 
be constant, already from f3 ~ 1.33, which is the 
transition point for the N T — 4. 

3. SUMMARY 

The phase diagram of the Bhanot-Cruetz || 
action in the fundamental and adjoint couplings 
has been crucial in understanding many proper- 
ties of the SU(2) lattice gauge theory, like the 
cross-over to the scaling region from the strong 
coupling region. Adding extra irrelevant terms to 
the action one obtains a modified action ([|) in 
which monopoles and vortices can be suppressed 
by setting the additional couplings to large values. 
Our finite size scaling analysis of L| -susceptibility 
for (3a — 0, Am = 1 and \e = 5, yielded 1.93 ± 
0.03 for the critical exponent u> = •y/v for lattices 
with N T = 4, thus verifying the naive universality. 
However, as a result of the suppression, the criti- 
cal coupling is shifted by about unity compared to 
the Wilson case. Our results on N T = 5 and 6 also 
yielded similar values for u> albeit with larger er- 



rors, confirming that the same physical transition 
was being studied this way as a function of the 
lattice cut-off, a. While the data for N^T 1 = aT c 
was found to vary slower with /? than expected 
from the asymptotic scaling relation for N T = 4- 
8, they did obey a similar relation with a factor of 
two larger exponent. A straightforward extrapo- 
lation suggests that the results from the modified 
action may merge with those of Wilson action for 
large N T ( of about ~ 25). Thus the suppression 
appears to make the transition from strong cou- 
pling to the scaling region much smoother than 
that for the unsuppressed Wilson action, allow- 
ing us to simulate the theory at smaller /3. It 
will be interesting to see if dimensionless ratios 
of physical quantities such as glueball masses or 
string tension with T c are constant in the range 
of critical couplings explored here. It will also be 
interesting to study such suppression in SU(3), 
and indeed SU(N) lattice gauge theories, since 
their phase diagrams are similar and the same 
mechanism is expected to work for them as well. 
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